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ABSTRACT 

We examine the cosmic growth of the red sequence in a cosmological hydrodynamic 
simulation that includes a heuristic prescription for quenching star formation that 
yields a realistic passive galaxy population today. In this prescription, halos dominated 
by hot gas are continually heated to prevent their coronae from fueling new star 
formation. Hot coronae primarily form in halos above ~ IO^^Mq, so that galaxies 
with stellar masses ~ IO^^'^Mq are the first to be quenched and move onto the 
red sequence at z > 2. The red sequence is concurrently populated at low masses 
by satellite galaxies in large halos that are starved of new fuel, resulting in a dip 
in passive galaxy number densities around ~ IO^^Mq that agrees qualitatively with 
observations. Stellar mass growth continues for galaxies even after joining the red 
sequence, primarily through minor mergers with a typical mass ratio ~ 20%. For the 
most massive systems, the size growth implied by the distribution of merger mass 
ratios is typically ~ 2x the corresponding mass growth, consistent with observations. 
This model reproduces mass-density and colour-density trends in the local universe, 
with essentially no evolution to z = 1, with the hint that such relations may be washed 
out by z ~ 2. Our simulation produces a high red galaxy fraction at both high galaxy 
overdensity, independent of stellar mass, and high mass, independent of overdensity, 
suggesting quenching mechanisms associated with both environment and mass; in our 
model, both are connected to the presence of surrounding hot gas. 
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1 INTRODUCTION 



A substantial fraction of galaxies today inhabit a tight locus 
in colour-magnitude space known as the red sequence. Such 
galajcies typically have elliptical morphologies, have little or 
no star formation, little cold gas, and live in dense environ- 
ments. T hey host the maj ority of stellar mass in the local 
universe (|Hogg et al.ll2002i i and therefore are a critical pop- 
ulation for understanding the global evolution of galaxies 
across cosmic time. 

Despite a growing wealth of data, the origin of red se- 
quence galaxies is still not fully understood. In the currently- 
favored cold dark matter (CDM) paradigm, massive galax- 
ies should be surrounded by halos of hot g as that can po- 
tenti ally cool and fuel new star formation (|White fc FrenkI 
Il99ll ). Indeed many red sequence galaxies are seen to have 
measurable X-ray emission, directly indicating hot, cooling 
gas. Yet these galaxies show little cold gas or ongoing star 
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formation. Debate continues as to what physical mechanisms 
can halt star-formation and keep it halted for much of cos- 
mic time with recent work favoring feedback from active 
galactic nuclei (AGN) as the primary energy source (e.g. 
ICroton et ai]|2006l : iHopkins et al.. 2008). But detailed ob- 
servations have revealed more subtle complications - for in- 
stance, a significant fraction of red-sequence galaxies have a 
"frosting" of current star fo rmation, though it contributes 
little to overall mass grow th (|Trager et al ]|l998l : lBower" et al.l 
1 19981 : pTrager et al]|2008l ). Furthermore, while most massive 
red sequence galaxies are ellipticals, pa ssive disks make up 
a large fraction at the low-mass end ijBundv et all |2010| : 
Ivan der Wei et alll2010h . Assembling all these observations 
into a coherent framework for passive galaxy evolution has 
proved challenging. 

More insight about the evolution of passive galaxies 
may be obtained by directly examining their history over 
cosmic timescales. This is now possible thanks to deep 
surveys that have placed constraints on the bright pas- 
sive galaxy populations out to 2: ~ 1 and even z > 
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2 (e . g. iBell et al.l 12004 iFaber et al.l l2007l: iBrown et alj 



20071; 


Cool et all 20081: Kriek et al. 


200a Stutz et al.l 20081: 


Tavlor ot al. 200f/: Brammer et al 


l200gl: iMarchesini et all 


2010: Whitakcr ct al. 2010). Ongoing surveys are ex- 



pected to probe passive galaxies with stellar masses 
down to ~ 1O^M0 (e .g. CANDELS; ICroein et all I2OIII : 
iKoekemoer et al] 120111 ). Current data indicates the bright 
end of the passive galaxy luminosity function has changed 
little (after correction for passive evolution) since z ^ 1, 
suggesting these galaxies were "in place" 7+ Gyrs ago. Fur- 
thermore, high-resolution imaging has revealed that passive 
galaxies are typically much smaller at z > 1, by factors 
^ 5 compared to their present-day de scendents (jPaddi et al.l 
l2005l : lvan Dokkum et ai1l2008l ■ l2010^ . These data likely pro- 
vide strong constraints on passive galaxy formation models, 
but such constraints have yet to be fully explored in the 
context of hierarchical models. 

Numerous physical processes are thought to contribute 
to the formation and evolution of passive galaxies, and many 
potential effects have proven difficult to rule out. Gas-rich 
galaxy mergers that induce a starburst and rapid black hole 
gro wth may drive the ini t ial quenching of mas s ive re d galax- 
ies (jSprineel et al.ll2005l : iHopkins et al.ll2006l . I2OO8I ). Alter- 
natively, hot gaseous coronae that form in halos above ~ 
10^^ Mq, supplemented by additional heating from an active 
galactic nucleus (AGN) or another energy source, may starve 
galaxies of fuel for star-formation (|Birnboim &: Dekell 20031: 
Keres et all l2005l: lOekel fc BirnboimI I2OO6I : ICroton et'all 
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tional possibilities may drive quenching: starvation and/or 
ram-pressure stripping by hot intra-group or intra-cluster 



actions with other galaxies (ICunn & GottI 19721: lAbadi et al. 


199S 


: Quilis et al' '2000*; 'Larson et al. 


Il980l: Bekki et al. 


2002 


:|Richstonc 1976: Moore ct al. 199^ 


. All these processes 



may act in combination to form passive galaxies as observed. 

Once "red and dead," the evolution of these galaxies 
may be deceptively complex. Although there is likely little 
or no new star formation, over its lifetime the stellar popula- 
tion sheds a significant fraction (~ 1/2) of its mass via stella r 
evolution ijjungwiert et al.ll200ll : iBruzual fc Chariot! l2003h . 
providing an additional local source of gas. Furthermore, 
passive galaxies are subject to mergers and accretion of small 
satellite galaxies, which alter their stellar populations and 
add mass. Mergers can also alter galaxy structure (e. g. from 
disky ordered rotation to spheroidal random orbit s) ([Barnes 
199C ), change the characteristic size of the galaxy lICox et al 



2006 ), strip stars int o the ICM (|Gallagher fc Ostrikeij| 19721 : 



physics of the transition to red and dead galaxies is less 
well constrained, and it is unclear which of the many physi- 
cal processes are driving this evolution. Cosmological hydro- 
dynamic simulations that dynamically track these complex 
physical interactions offer a maturing framework for eluci- 
dating the key physical drivers, and simulations have now 
begun to explore the evolution of the passive galaxy pop- 
ulation, particularly in connection with blac k hole growth 
llDi Matteo et al.1 l2005l : iHopkins et al] I2OO6I : ISiiacki et all 
l2007h . 

In this work, we study passive galaxy growth from 
high redshift until today using cosmological hydrodynamic 
simulations. These simulations incorporate physically- 
motivated but heu ristic quenching mechanisms (described in 
iGabor et aLll201ll ) that yield a 2 = population of passive 
galaxies whose colour and luminosity distributions match 
observations. In particular, we include a prescription where 
we ensure that halos dominated by hot gas are continuously 
heated to keep circum-galactic gas hot. This simple and 
extreme prescription, approximating the effects of a heat 
source such as an AGN radio jet, successfully cuts off the fuel 
supply for star formation in massive halos whi ch eventually 
yields passive galaxies. In contrast, we found in lGabor et alJ 
that a quenching prescription based on major merg- 
ers alone fails to produce the observed properties of pas- 
sive galaxies. Although significant challenges remain for our 
favoured model, it provides a general qualitative guide for 
models where hot massive halo s are the mai n drive rs of pas- 
sive galaxy formation. While iGabor et al.l l|201ll) focused 
on the z = Q red galaxy population and compared various 
quenching mechanisms, here we focus on our most successful 
quenching mechanism and its implications for passive galaxy 
evolution from high redshift until today, including environ- 
mental factors. Our model should not be considered physi- 
cally correct in detail, but rather as broadly illustrative of 
the impact of some (unspecified) quenching mechanism that 
keeps hot halo gas hot on passive galaxy evolution over the 
past 10 billion years. 

We begin by presenting an overview of our simulations 
in 321 Then we present our results, beginning in §3] where 
we examine number density and colour evolution, §|3] where 
we look at how passive galaxies grow in mass and size via 
mergers, and §[S] where we study mass and colour evolution 
versus environment. Finally, we summarize and discuss our 
conclusions in § |6l 



Murante et alll2004l ). and possibly add cold gas. Such pro- 2 SIMULATIONS 



cesses may all contribute to the observed evolution of the 
red sequence, and hypotheses addressing the observed lack 
of evolution must explain why these processes did not act 
or conspired to yield little evolution. 

From a theoretical point of view, the challenge of galaxy 
formation is to assemble a collection of physical processes 
within a hierarchical structure formation scenario that yields 
the observed properties of both star-forming and passive 
galaxies at all redshifts. Recent work on star-forming galax- 
ies suggests that their growth can be well-described as being 
fed by cold, filamenta ry accretion and r egulated by strong 
galactic outfiows (e.g. iDave et~al ] |2011bl lah. Such processes 
are either observed and/or physically well-motivated. The 



2.1 Simulation methodology 

We analyze the same co smological hydrodynamic simula- 
tions described in lGabor et al. (2011). For completeness, we 
summarise the simulations here. They were run with an ex- 
tended version of the N-b ody -I- smoothe d particle hydrody- 
namics code GADGET-2 (|Springelll2005l '). In addition to the 
basic N-body and hydrodynamics calculations, our version 
of the code includes sub-resolution modelling of gas cooling, 
star-formation, a model for chemical enrichment via AGB 
stars and both Type la and core-collapse supernovae, galac- 
tic winds associated with star-formation, and simple pre- 
scriptions for quenching star-formation. 
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We include both primordial and metal-line 
cooling assuming col lision al ionisation equilibrium 



(jSutherland fc Dopital \l99± . with the metallici- 
ties self-consistently tracked wit hin simulations; see 
lOppenheimer fc Pavel l|2006l . |2008| ) for details. We in- 
clude h eating from a metaaa lactic photo-ionising back- 
ground ijHaardt fc MadaullioOll '). assuming all particles are 
optically thin. We truncate the cooling curve at 10* K be- 
cause our simulations lack the resolution to properly model 
the interstellar medium; below this temperature, particles 
can only be heated or cooled via adiabatic expansion. 

For star formation, we employ the two- phase model of 
ISpringel fc Herngu ist (2003) based on the analytic descrip- 
tion of lMcKee fcOstriker (,19771 ). Gas particles above a den- 
sity threshold of 0.13 cm~^ are treated as cold, star-forming 
clouds embedded within a hot diffuse medium. The masses 
in the cold and hot phases are determined by the balance 
between evaporation by core-collapse supernovae and con- 
densation. The star formation rate is then determined by the 
timescale over which cold clouds collapse into stars, and the 
model parameters are set so that the resulting star forma- 
tion rates ar e consistent with the observed relation to surface 
gas density (|Kennicuttlll998l ) . Star- forming gas particles are 
converted into collisionless star particles stochastically, with 
a probability derived from the star formation rate. Gas par- 
ticles may spawn up to two such star particles, conserving 
total mass. 

As a gas particle is undergoing star formation, it self- 
enriches with metals from core-collapse supernovae. Fur- 
thermore, star particles share energy, mass, and metals 
with neighboring gas particles as a result of stellar mass 
loss from AGB stars and Type la supernova e. Stell ar mass 
loss is calculated by assuming a IChabrieJ (|2003l ) initial 
mass function, applying the stellar population models of 
iBruzual fc CharlotI (|2003l . , BC03), which account the mass 
lost by a stellar population at discrete times after the ini- 
tial s tar- formation event . We u se Type la supernova rates 
from IScannapieco et al] (|2006l ). and each such supernova 
results in the production of metals (mainly iron) that are 
shared with neighboring gas particles. We note that in simi- 
lar simulations, the predicted metallicities in hot gas around 
mas sive galaxies are in good agreement with observed val- 
ues (|Dave et al.ll2008l ). 

Our code further includes feedba ck in the form of galac- 
tic w i nds d riven by star-formation (jOppenheimer fc Pavel 
l2006l . l2008t ). Just as a star-forming gas particle has some 
probability of being converted into a star particle, it has a 
probability of being kicked in a wind, which is given by 77 (the 
mass loading factor) times the star formation probability. A 
wind particle is expelled from its host galaxy at a velocity 
Vw typically a few hundred km s~^, and it is briefly decou- 
pled from normal hydrodynamic interactions to approximate 
chimneys of escaping gas (Springcl & Hcrnquist 20d^. The 
velociti es are chosen t o match observatio ns of local galaxy 
winds (jMartinl l2005l : iRupke et al] |2005|) and also match 
those seen at higher redshifts (e.g. Steidel et al.ll201ol ). The 
mass loading factors scales as the inverse of the galaxy 
circular veloc ity, which is calculated fr om an on-the-fly 
galaxy finder (|Oppenheimer fc Davdl2008l ). and follows scal- 
ings a s predicted for momentum-driven winds (jMurrav et al.l 
I2OO5I ). With this prescription, our simulations match a broad 
array of observational constraints on star-forming galax- 



ies and the intergalactic medium (lOppenheim er fc Pavel 
2006l.l2008l:rOppenheimer et al.ll201Gl:|Pave et al.l l2006. 2o6i; 



Finlator et al]|2007l : iFinlator fc Pavell2008l ). 

Pespite their broad success, winds driven by star- 
formation generally do not result in the formation of red 
and dead galaxies. We incorporate an additional heuris- 
tic quenching mo del specifically to solve this problem 
(Gab or et al.ll201ll ). For this model, we run a spherical over- 
density algorithm on-the-fly to identify galaxy halos and sep- 
arate gas above and below 250,000 Kelvin in each galaxy's 
halo. In halos with > 60% of all gas above this tempera- 
ture cutoff, we apply constant heating (at every time-step) 
to all circum-galactic halo gas to force it to remain at the 
halo virial temperature. Such heating could plausibly re- 
sult from an AGN radio jet, but our model in fact uses 
more energy than thought to b e available from observed 
AGN sources (Gabor et al. 201^. We exclude star-forming 
gas particles from this heating since a radio jet is unlikely 
to disrupt dense star-forming clouds on galactic scales. The 
heating we apply to circum-galactic gas is sufficient to pre- 
vent it from condensing onto the galaxy, thus starving the 
galaxy of new fuel for star-formation. This model results in 
a bimodal colour distribution of galaxies, and red galaxy 
luminosity functions that match observations of the local 
universe. Although the model has some difficulties (like the 
required energy mentioned above), it should be a good rep- 
resentation of models where galaxies are quenched due to 
starvation enabled by their surrounding hot coronae. 



2.2 Runs and galaxy identification 

For this paper, we employ a simulation of a 48h^^ comoving 
Mpc random cosmological cube with 256^ dark and 256'^ 
gas particles that incorporates all the above physics. We 
use a Wilki nson Microwave Ani sotropy Probe concordance 
cosmology (jKomatsu et al.ll2009^ with Hq = lOO/i = 70km 
s~^ Mpc~^, matter density Q.m = 0.28, baryon density — 
0.046, a cosmological constant with Q,a ~ 0.72, root mean 
square mass fluctuation at separations of 8 Mpc erg — 0.82, 
and a spectral index of n = 0.96. Our simulation uses a 
gravitational softening length of 3.75/i~^ kpc. The initial gas 
particle mass is 1.2 x 10* Mq, the typical star particle mass 
is half that, and our simulation results in ~ 3000 resolved 
galaxies at z = 0. 

We save snapshots of the simulation at 108 redshifts, 
starting at 2; = 30 and ending at z = 0. The time between 
snapshots ranges from a few tens of Myrs at high redshift to 
~ 300 Myr at low redshift. The snapshots contain informa- 
tion for simulation particles, such as position, velocity, mass, 
metallicity, gas density, gas temperature, and star-formation 
rate. From these particle data, we determine galaxy proper- 
ties to compare with observables. 

We use Ski d El t o identify galax ies (cf. 
iGelb fc BertschingeJ Il994l : iKeres et all I2OO5I ). Skid 
provides a list of member particles (star and star-forming 
gas) for each simulated galaxy. The sum of member star 
particle masses is then the galaxy stellar mass, and the 
instantaneous star formation rates of the gas particles are 
summed to give the star formation rate of the galaxy. 
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We then calcu l ate ga laxy s pectra using th e mod els of 
iBruzual fc CharlotI (|2003l ). as in iFinlator et al] l|2006l ). We 
treat each star particle as a single stellar population with an 
age and metallicity determined directly in the simulation. By 
adding up the spectra of all star particles within a galaxy, we 
obtain the spectrum of that galaxy, from which we measure 
galaxy colours and magnitudes in various bands. Though we 
focus on the passive galaxy population that is thought to 
be mostly dust-free, when considering star-forming galaxies 
we some times include a presc ription for dust obscuration 
following iFinlator et alj l|200el '). 

2.3 Building merger trees 

Beyond knowing galaxy properties at a given redshift, we 
wish to study the histories of individual galaxies as they 
evolve. For this we must connect each galaxy at z = to its 
progenitor galaxies at earlier redshifts: i.e., build a merger 
tree. We do so by determining, for every star particle in a 
given galaxy, which galaxy from an earlier snapshot that star 
particle lived in. If the star formed recently it wouldn't have 
lived in any previous galaxy. By tracing star particles in this 
way, we determine the immediately preceding progenitors of 
each galaxy. Galaxies with multiple progenitors must have 
undergone a merger (or at least some galaxy interaction, 
such as stellar stripping). 

Creation of a merger tree from a cosmological simu- 
lation involves handling a number of pathological circum- 
stances owing to the stochastic nature of the gas and star 
particles. A common annoyance is that if two galaxies fly by 
each other in a close encounter. Skid (and most other galaxy 
finders) will often identify them as a single galaxy during one 
or a few timesteps, but then again as separate galaxies at 
later times. More rarely, a group of stars that make up a 
galaxy at one timestep may be assigned to no galaxy at all 
in a subsequent timestep (i.e. they are part of the intraclus- 
ter light), giving the appearance that a galaxy disappeared. 
This can happen with quite large galaxies (> 1000 particles). 
While these issues do not affect the majority of galax;ies and 
hence do not strongly impact the overall results, they must 
nevertheless be handled in an appropriate manner. 

To su rmount these difficu lties, we follow strategies de- 
scribed in iMaUer et all (j2006l l. who dealt with the same is- 
sues. We distinguish Skid groups from galaxies, the latter 
being defined in terms of our merger trees. For each star 
particle in a Skid group at 2 = 0, we trace the history of its 
host Skid groups over all timesteps. If the star particle does 
not appear in any Skid group during some timestep, then 
we assign it to a virtual group where it is the only member. 
We then assign each group (real and virtual) to a galaxy. At 
the highest-redshift timestep of interest, each Skid group is 
defined as a separate galaxy. In the subsequent timesteps, 
each group is assigned to the same galaxy as its most massive 
progenitor from the immediately preceding timestep. Thus 
if two separate Skid groups share the same progenitor, then 
they are both assigned to the same galaxy, even though Skid 
would identify them as separate. This helps handle the case 
where two merging galaxies undergo a close passage, then 
separate before the final merger. Note that the precise tim- 
ing of final coalescence may be incorrect by up to ~ 1 Gyr, 
but our results do not rely on this timing aspect. Likewise, 
virtual groups that were not part of any Skid groups are as- 



signed to galaxies that they once inhabited. Once all groups 
are assigned to galaxies, their constituent particles can also 
be assigned to those galaxies. 

The above prescription assumes that once galaxies are 
joined by Skid, they will eventually merge. But some fiy-by 
interactions will not lead to mergers. To account for this, 
we identify galaxies at 2; = that are composed of multiple 
Skid groups and separate those groups into distinct galax- 
ies. We find the earliest timestep when both groups' stars 
are assigned to the same galaxy, and from then on assign 
a new galaxy ID to those stars in the smaller of the two 
groups. With every star particle assigned to a galaxy at ev- 
ery timestep, we then identify all progenitor galaxies to each 
galaxy. 

2.4 Environment measures 

In order to study where red galaxies emerge in our sim- 
ulation, we measure the local galaxy density around Skid 
galaxies. We use two density measures. An intuitive and sim- 
ple approach is to count the number of si mulated galaxie s 
within lh~^ Mpc of the galaxy of interest (|Blantonll2006l ). 
For some purposes this measure too noisy, as it places galax- 
ies into discrete bins of density, and some galaxies have few 
neighbors within Ih^^ Mpc. Thus, as an alternative we use 
Smooth[f|, which calculates the mean density at each galaxy's 
location after smoothing over the nearest 16 galaxies with a 
spline kernel (similar to an SPH density calculation). This 
allows us to quantify low densities less stochastically than 
by simply counting galaxies within some fixed distance. In 
some cases we use a galaxy's local overdensity to compare 
with observations: the overdensity is measured relative to 
the mean resolved galaxy density in our simulation volume. 



3 EVOLUTION IN COLOUR AND NUMBER 
DENSITY 

3.1 Color- mass diagrams and stellar mass 
functions 

Recent observations enable two simple but powerful con- 
straints on the emergence and evolution of passive galax- 
ies: their colours and number densities. In this section we 
examine how colours and number densities (parameterized 
via the galaxy stellar mass function) evolve with redshift in 
our simulations. We note that our simulations do not per- 
fectly match present-day observa tions of passive galaxies, 
as detailed in iGabor et all l|201lD : for instance, the colors 
are slightly too blue and the red sequence does not have 
the correct slope (likely a metallicity effect). Nonetheless, 
our simulations clearly identify two general trends in pas- 
sive galaxy evolution: (1) the first passive galaxies emerging 
at 2 > 2 are at the massive end of the galaxy stellar mass 
function, followed shortly by small satellites, then finally the 
middle region of IO^'^A/q galaxies. (2) The most massive red 
galaxies in our model grow substantially at late times. 

Figure [T] shows colour-mass diagrams and galaxy stellar 
mass functions for simulated galaxies at redshifts 0, 0.3, 1, 
and 2. The g—r colours are rest-frame SDSS filter bands, and 

^ http:/ /www-hpcc. astro. washington.edu/tools/smooth. html 
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Figure 1. Rcdshift evolution of colour- mass diagrams (top row) and galaxy stellar mass functions (bottom row). Redshift increases 
from left to right, as labeled. Red and blue galaxies are separated using a redshift-dependent cut, as illustrated by the red and blue 
colour-coding in the CMD. In the mass function panels, solid lines denote the simulated red sequence, dashed lines the simulated blue 
cloud, and observed data points are taken from the literature. The hatched region indicates poorly resolved galaxies. Massive red galaxies 
are present at 2 = 2, but not in sufficient numbers compared to observations. There is a dearth of intermediate mass galaxies (~ 10^° Mq) 
that is more pronounced at higher redshifts. 



do not include any correction for obscuration due to dust - 
these are intrinsic colours. We divide the galaxy population 
into blue and red using a straight line i n q — r colour versus r 
magnitude space in the same way as in lGabor et al.l l|201ll ). 
This translates to a similar separation in colour-mass space, 
as indicated by the sharp demarcation between blue and 
red points in the top left panel. We evolve the normaliza- 
tion (but not the slope) of our colour separation in a simple 
linear way that approximately accounts for passive evolu- 
tion: ysepiz) = J/sep,z=o — 0.1 * z. Here j/sop is the y-intercept 
(i.e. colour) of the line of separation, and z is the redshift. 
The resulting separation indicated by the colour-coding of 
blue and red galaxies in the top panels of Figure [1] suggests 
this is a reasonable choice. The hatched region on the left of 
the figure denotes the mass range of galaxies that may be 
inadequately resolved. 

The colour evolution of the galaxy population is clear 
from the figure. The main feature is the strong growth in 
the passive galaxy population from z = 2 — !> 0. The blue 
galaxy population also becomes redder with time, at roughly 
the same rate as the red sequence. This indicates that even 
galaxies in the blue cloud have their colour evolution driven 
by passive stellar aging. Recall that we do not include ex- 
tinction here, but particularly at high redshifts, many of the 
blue galaxies may lie in the red region had we included dust; 
this complicates comparisons with observations, which typ- 
ically have difficulty separating these populations (passive 
vs. dusty star-forming) at higher redshifts. 

The lower set of panels of Figure [1] shows galaxy stel- 
lar mass functions separated into blue (dashed line) and 
red (solid line) pop ulations. We plot mass function obser - 
vati ons taken from iBeU et al.1 (|2003l ). iDrorv et all (|2009t ). 
and iBrammer et al.1 (120111 )" who likewise split galaxies into 
star-forming and quiescent populations (red vs. blue points), 
though with slightly different criteria. We note that our sim- 



ulated mass function has too few massive galaxies, but ap- 
proximately the correct shape at z = (increasing all stellar 
masses by ~ 50% provides a good match). Although this 
model provides a goo d match to the lum inosity function of 
red galaxies at 2 = (iGabor et al.ll201ll ). the corresponding 
masses of these galaxies are too small, indicating a problem 
with their mass-to-light ratios. Increasing the hot gas fracion 
above which we apply quenching from the value of 60% used 
here can effectively shift the mass function to higher masses 
(|Gabor et al.ll201Gl ) and may provide better agreement with 
this data, though likely at the expense of the luminosity 
function agreement. These small differences notwithstand- 
ing, we are mostly interested in the evolution of the red 
sequence, and hence we focus on qualitative aspects of this 
model rather than tuning to match observables. 



From the figure we see that the simulated blue galaxy 
stellar mass function evolves only weakly with redshift, but 
the red galaxy mass function obviously undergoes drastic 
growth from 2 > 2 to 0. The rapid buildup of passive galax- 
ies refiects the quenching o f star-forming galax ies due to 
starvation in our model (see iGabor et alI(|201ll )V The red 
sequence grows from mostly non-existent at 2; ~ 3 to domi- 
nating the galaxy population at 2 ~ 0. The mass functions 
clearly show that the most massive galaxies are the first to 
become red, followed by low-mass galaxies that turn out to 
be satellites to the massive ones. This directly results from 
the fact that the most massive galaxies in our simulations 
are the first to form a corona of hot gas. Massive galaxies 
tend to live in massive dark matter halos, and massive dark 
matter halos have low enough densities that virial shock 
heating outpace s radiative cooling, allow ing a hot gaseous 
halo to develop l|Birnboim fc Dekei|[2003l ) . We note that our 
simulation volume is limited, so that we are likely under- 
estimating the prevalence of the most massive systems that 
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would be the very first ones to become passive at the highest 
redshifts. 



3.2 A dearth of passive galaxies at M* 



10^°Mq 



Examining the red sequence growth in more detail, we find 
that the first passive galaxies form at the massive end by 
~ 3, which is likely a lower limit since our limited simu- 
lation volume does not capture the most massive structures 
in the Universe. Soon thereafter, small satellites of the first 
quenched galajcies become red as well, owing to starvation of 
gas in growing hot halos (jSimha et al ] |2009h . Intermediate- 
mass galaxies with Mateiiar ~ IO^^Mq are the last to quench, 
filling in the gap between the centrals and satellites such 
that by 2 ~ this bimodality within the red sequence is not 
so evident. Therefore, a strong prediction of this quenching 
mechanism is that there is a distinct gap between the mas- 
sive and low-mass ends of the red sequence that becomes 
more pronounced with redshift. 

Galaxies around IO^^A^q at z > 1 are typically nei- 
ther massive enough to support their own hot coronae, nor 
close enough to larger galaxies to live in the hot gas en- 
vironment of a (proto-)group or cluster. By redshift zero, 
however, many galaxies within this mass range live in halos 
dominated by hot gas due to the presence of very massive 
neighbors. As we will discuss further in § [S] at times they 
may not even occupy the same virialised halo, but rather 
live in a nascent super-structure with pervasive hot gas that 
triggers quenching. 

The increasingly prominent dip in the red galaxy mass 
function at M* ~ IO^^Mq to higher re dshifts is consistent 
with existing observational constraints. IStutz et al.l (|2008| ) 
found six massive (> IO^^'^Mq) but zero low-mass galax- 
ies with old stellar populat ions at 2 ^ 3 in t he Hub- 
ble Ultra-Deep Field fHUDF: [Beckwith et allbood ). despite 
being sensitive t o stellar masse s below 10^° Mq ( see also 
iToft et al.|[2005h . iRudnick etHI (2009') noted that in clus- 
ters, the low mass end of the mass function must build 
up more rapidly than the massive end from z ^ 0.8 to 
2 = in ord e r to " catch up" in number density. In the field, 
iDrorv et all (|2009l ) traced the evolution of the red galaxy 
stellar mass function for z < 1, finding a dip at medium 
masses that is more prominent at higher redshift; some of 
these data are plotted in F igure [l] Similar tren ds have been 
noted at low redshift (e.g. iBaldrv et al]|2008l ). Our simple 
model, where hot gas dictates the formation of the red se- 
quence, provides a compelling physical picture behind this 
trend, which we describe next. 



3.3 Central and satellite evolution along the red 
sequence 

The separation in mass range between central and satellite 
galaxies along the red sequence is more pronounced at ear- 
lier epochs, leading to a sort of bimodality along the red 
sequence itself. This is naturally understood within a hier- 
archical structure formation model, as we outline here. 

In our model, the initial emergence of the massive 
red sequence is directly tied to the requirement that a 
galaxy's halo be dominated by hot gas for quenching 
to begin. Early spherically-symmetric analytic models of 
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Figure 2. Top: Rate of change of red galaxy number density 
with respect to redshift, for three stellar mass bins (solid, dotted, 
and dashed l ines). We compare to smoothed derivatives based 
on data from iBrammer et al. I ll201ll , , squares and diamonds for 
two mass bins). At z < 0.5 in our model, the number of low- 
mass galaxies starts to catch up to the number of more massive 
(middle bin) galaxies, as the derivative there is higher for low- 
mass galaxies. Bottom: Red galaxy fraction as a function of z 
for the same three bi ns in stellar mass, again with data from 
iBrammer et al. I l l201ll 'l. Our model under-predicts the observed 
red galaxy fraction at Al < lO"'^'''^. 



galaxy evolution assumed that gas in the halo of a forming 
galaxy should be vir ialised ()Rees fc Ostrikeij|l977l : ISilklll977l : 
IWhite fc Re5ll978l i. It was realised, however, that gas in 
low-mass halos will have short cooling times, s o that it never 
reaches the virial temperature of the halo jBinnevI Il977l : 
I White fc Freni3ll99ll : iBirnboim fc Dekelll2003l i. Several re- 
cent hydrodynamic studies have shown that a hot circum- 
galactic corona only f orms, via the virial shock, in halos 
abov e a few xlO^M© ()Keres et al.ll2005l : IBirnboim fc Dekel 
l2003h . These analyses neglected metal-line cooling, which 
is particularly important in massive hot halos that are 
heavily enriched (Dave et al.ii200S ). Our simulations includ- 
ing metal-line cooling indicate a hi gher mass threshol d of 
~ IO^^Mq for a hot corona to form (jCabor et al.ll2010l '). In 
accord with previous analyses, we find that this threshold 
does not evolve significantly with redshift. 
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The key assumption in our model is that in such hot- 
gas dominated halos, gas is very inefficiently deposited 
onto the galaxy, likely owing to the contribution of some 
preventive feedback process(es) like AGN heating (e.g. 
Dekel fc Birnboirall2006l : ICroton et al.ll2006l: ICattaneo et all 



20061 : iBower et al.l[200alBirnboim et al.ll2007C Since galaxy 
stellar masses are well-correlated with halo masses below the 
cutoff for quenching, only central galaxies with the highest 
stellar masses will live in hot halos at early times. 

As the universe evolves, more galaxy halos become large 
enough to support hot coronae. Massive central galaxies thus 
continue to move onto the red sequence. Our hot gas frac- 
tion cutoff of 60% roughly corresponds to a stellar mass of 
about 10^°'*Mq, though there is significant scatter. Once 
galaxies exceed this stellar mass, they will tend to migrate 
toward the red sequence over the span of 1 — 2 Gyr. We 
note that this mass cutoff is unlikely to be correct at 2 > 2, 
since there are seen to be a number of star-forming galax- 
ies seen with stellar masses app roaching and ev en exceeding 
10"M fl (e.g. lDaddi et"lLll2007l ) . Simulations bv lDekel et al.l 
l|2009[) suggest that at these redshifts, dense cold inflow can 
penetrate hot halos and therefore fuel star formation even 
in the presence of a hot corona. At lower redshi fts, this phe- 
nome non is expected to weaken considerably l|Keres et al.l 
I2OO9I ). and thus our quenching model may be more appro- 
priate for 2 < 2 or perhaps z < 1.5 passive galaxies. 

What happens with satellite galaxies? T he halo occu- 
panc y of small halos is typically unity (e.g. iBerlind et al.l 
I2OO3I ). that is, small halos tend to be dominated by a sin- 
gle central galaxy with only very small satellites present. 
This is because the cold accretion that feeds galaxy growth 
tends to bypass satellites, thereby disproportionately grow- 
ing the central object. When such a central galaxy becomes 
quenched, it satellites also undergo quenching by starvation 
owing to the growth of the same hot halo that triggers our 
quenching mechanism. But these satellites are quite small, 
and therefore initially populate the very low-mass end of the 
red sequence, leaving a strongly bimodal mass distribution 
of passive galaxies. 

These quenched halos then merge with nearby halos 
owing to hierarchical growth. As the smaller halo falls in, 
its central galaxy becomes a satellite and is quenched b y 
the hot ICM on a 1 - 2 Gyr timescale (jSimha et al.ll2009l ). 
These infalling (former) central galaxies are typically larger 
than the original satellites within the quenched halo, and 
hence begin to populate the intermediate-mass portion of 
the red sequence. They of course also bring in additional low- 
mass satellites. Therefore as time progresses, hierarchical 
growth combined with starvation by the hot corona results in 
a gradual "filling in" of the intermediate-mass red sequence. 
These trends will be explored in more detail, particularly in 
relation to environment, in ^ 

Recent semi-analytic models bv iDe Lucia et all (I2OIII ) 
analogously found a bimodality in elliptical galaxy masses 
with a minimum at Mgteiiar ~ IO'^'^Mq. They argue that this 
arises in their model owing to mergers probabilities increas- 
ing with mass and disk instabilities contributing to bulge 
growth. Our simulations yield qualitatively similar results, 
but in our case passive galaxy formation is not driven by 
mergers or instabilities, but rather by starvation by a hot 
halo. It is worth noting that their model is geared towards 
producing elliptical morphologies, while ours is focused on 



producing passive systems, which are highly overlapping 
population in the real Universe but may arise from different 
phys ical mechanisms (and not exactly coevallv: lBundv et al.l 
I2OIOI ). It would be interesting to compare and develop dis- 
criminating predictions for these scenarios; we leave this for 
future work. 



3.4 Red galaxies undergo strong evolution 

Figure[2]illustrates the evolutionary trends in passive galaxy 
number density more explicitly. We consider three bins in 
stellar mass, and show how the number density of passive 
galaxies in each bin changes with redshift via the derivative 
drired/dz. Passive galaxies with masses between 10'^° '' and 
10^° '^ Mq are the most numerous at all epochs, and expe- 
rience the highest growth rate over most of the history of 

0' 



the universe. Less massive galaxies with Afstoiiar ~ 10 
i.e. within the "gap" in the mass function, emerge later but 
their growth rate accelerates to later epochs, being higher 
than more massive galaxies at z < 0.5. This illustrates the 
catch-up growth in the numbers of middle mass galaxies. 
Finally, red sequence galaxies more massive than 10" M0 
are formed owing to migration from the lower-mass bins as 
they acquire mass through mergers, which is a relatively 
slow process, and hence their growth rate is the slowest. 

Comparison betw een simulated a nd observed red 
galaxy mass functions l|Bell et all l2003l : iDrorv et al.l l2009l : 
iBrammer et al.l [loij) in Figure [T] shows that our model 
under-predicts the numbers of massive (> IO^^Mq) red 
galaxies, increasingly so at higher redshifts. In the real Uni- 
verse, significant populations of massive red galaxies exist at 
z > 2, and they grow appreciably in number sin ce that epoch 
(but growth is modest from z ~ 1 to ~ 0.3 Brown et al 
'20071: ICool et al.1 l2008l : iTavlor et al] l2009l : iBrammer eFaT 
20111 ). To quantify the discrepancy between observations 
and our simulation in Figure [5J we also show smoothed ob- 
servational number density growth rates for two mass bins 
derived from th e abso lute number density measurements of 
IBrammer et al.l l|201ll . squares and diamonds) . This confirms 
that the growth rate of the largest passive galaxies is some- 
what too low in our simulations. 

The bottom panel of Figure [5] shows the fraction of 
galaxies that are red in different mass bi ns as a function 
of red shift. Again we include data from IBrammer et al.l 
( l201lf ). At z < 1.5 our simulation matches the observations 
well for the most massive galaxies only. The simulation un- 
derpredicts the red fraction for galaxies less massive than 
10^" '^ Mq, probably due to an excess of blue star-forming 
galaxies at these lower masses. 

To summarise the last few sections, our model pre- 
dicts that galaxies with Afstoiiar ~ 5 x IO^^Mq are the 
first to become passive, but it underpredicts the numbers 
of galaxies with Mstoiiar > 10^^ Mq. Galaxies in the "gap" 
at Afgtciiar ~ lO^'^AfQ catch up with time. The key feature 
driving these trends is that our quenching mechanism, tied 
to the hot gas coronae that form in halos ^ lO'^^AfQ, selects 
a characteristic mass for passive galaxies. That is, once a 
star-forming galaxy achieves a mass > 1O^°-^M0, it wiU be- 
come quenched and move to the red sequence. Less massive 
galaxies are only quenched as satellites (or nearby halos) 
of more massive halos. More massive galaxies can only be 
obtained via merging. While quantitative details may vary 
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Figure 3. Example evolutionary paths in the colour-mass dia- 
gram. We show the paths of three galaxies (coloured lines) as they 
evolve through the colour-mass diagram, as represented by z = 
galaxies (small points) in the background. T his is analogous to 
the schematic diagram of iFaber et al. 1 1 I2OO7I) . Both background 
galaxies and paths include a prescription for dust obscuration, 
and paths are corrected for passive colour evolution. Black lines 
indicate redshifts along the evolutionary paths, as labelled for the 
most massive galaxy. Galaxies grow in stellar mass along the blue 
sequence, then move vertically to the red after quenching. Once 
on the red sequence, galaxies grow only via mergers and accretion 
of satellites. Such growth can be substantial, as for the path of 
the most massive galaxy shown. 

with parameter choices, these are generic outcomes for hier- 
archical quenching models keyed to a critical halo mass. 



3.5 Paths to the red sequence 

Here we focus on the path of individual galaxies in colour- 
mass space: Given a red galaxy at z ~ 0, what path did 
it take to get there? It turns out that massive galaxies in 
our models follow path s muc h like those in the schematic 
diagram of iFaber et al.l (|2007| ): they move onto the red se- 
quence at masses near the turnover in the mass function, 
then grow further through mergers and accretion of satellite 
galaxies. 

In Figure |3] we show the paths of three typical galax- 
ies (coloured lines) in the colour-mass diagram. We plot a 
background of simulated galaxies at z = (points), and we 
correct for passive evolution of colours for the galaxy paths 
in the same way that we evolved our line separating blue 
and red galaxies H3.ll Here we also include a simple pre- 
scription for dust that is tied to the star-formation rate of 



each galaxy ( 


Somerville et al.ll200ll:IWanE & Heckmanll996l 


iGabor et al.l 


2010h. 



All three of these galaxies build stellar mass via star- 
formation for several Gyrs before moving onto the red se- 
quence. Once quenched, they stop building stellar mass and 
cross vertically and rapidly over the green valley to the red 
sequence in < 2 Gyr (see the star-formation histories in 



iGabor et al]|201ll ). The two less massive galaxies make this 
transition at late times (2 < 0.5), and do not change in 
mass or colour after reaching the red sequence because they 
do not have time to undergo mergers. 

The most massive galaxy, in contrast, reaches the red 
sequence at z ~ 2 when its mass is a few xIO'^'^Mq, and 
then continues to grow in mass along the red sequence via 
mergers, eventually obtaining a mass of ~ 10" M0. The 
colour evolution along the red sequence is generally towards 
the blue (when corrected for stellar evolution), by virtue 
of accretion of smaller galaxies that are typically less red 
owing to lower metallicities. This is the basic reason why, 
as we showed in iGabor et all l|201ll ). the predicted red se- 
quence slope is too shallow and does not match the observed 
red sequence slope, particularly at the massive end. Given 
that our star-forming galaxies li e on the mass-m etallicity 
relation IFinlator & Dave 2008; Dave et aLlbOllal ). it is dif- 
ficult to envision ways to avoid this natural consequence of 
merging. One possibility is that the amount of merging is 
overestimated in our models, perhaps because interactions 
in the real Universe result in satellites being stripped into 
the ICM rather than joining the central galaxy. If so, intr- 
acluster light is predicted to have significantly lower metal- 
licities than the central galaxy; this is difficult to test obser- 
vationally at present. Another possibility is that there are 
too many low-mass galaxies at early epochs in our mod- 
els, as c omparisons to high -z mass functions preliminarily 
suggest l|Dave et al.ll2011bl ). which provides excess fodder 
for late-time merging. Clearly the detailed evolution of the 
most massive systems is not fully vetted in our current sim- 
ulations. 

More broadly, these paths are fairly representative for 
galaxies in their respective mass ranges. Low-mass galaxies 
tend to move onto the red sequence at late times, and grow 
little after doing so, while massive galaxies become passive 
at earlier epochs and grow in mass by merging with smaller 
galaxies. 

3.6 Specific star formation rates 

A key observational constraint on the transition of galax- 
ies onto the red sequence is provided by the distribution 
of specific star formation rates (sSFR) versus stellar mass. 
This contains similar information to the colour-mass plots 
above, but here we additionally compare our current quench- 
ing simul ations with previou s similar runs without quench- 
ing from iDave et al.l (|2011lJ ). We particularly examine the 
stellar mass range over which the transition to passive oc- 
curs, since the observed lack of a sharp delineation in stellar 
mass between passive and star-forming galaxies has been 
used as an argument against a quenching prescription based 
on halo mass (as ours is, approximately). 

Figure |4] shows sSFR as a function of stellar mass for 
a simulation without any quenching (green points) and the 
halo quenching run presented in this paper (red) . Blue lines 
show a running median for each. Galaxies with sSFR^J 
10"'' Gyr~^ are shown along bottom of the plot. The no- 
quenching run employs the same cosmology, volume, input 
physics (besides quenching), and outflow model. But it has 
more particles (2 x 384^) , and consequently has somewhat 
better mass and spatial resolution, and hence its resolved 
galaxy population extends to slightly lower masses. 
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Figure 4. Specific star formation rate (sSFR) as a function of stellar mass at redsliifts (left), 1 (middle), and 2 (right). Red points 
correspond to our halo quenching model and green points to a simulation without quenching. Blue solid lines show the running median 
sSFR for both models. Many quenched galaxies have no star-formation, so we assign a minimum sSFR of 10-3 Gyr-l. Quenching has 
little impact on the sSFR-Af, relation below ~ IO^^'^Mq, but causes a drop in sSFR above this threshold at all redshifts. The transition 
quenching mass is blurred by satellite galaxies ~ IO^^Mq becoming quenched. 



The form and evolution of the sSFR fo r the non- 
quench ed case is discussed extensively in iDave et al.l 
l|2011bl ). Key points are that the relation is fiat at hieh-^ 
moving towards a mildly negative slope by z ~ 0, sSFR in- 
creases with redshift at a given stellar mass steadily out to 
2 = 2, and the scatter is typically ~ 0.2 dex or less. The 
relation fundamentally arises from gravitationally-driven 
smooth cold accretion that dominates galaxy fueling , and 
diminishes with time as the Universe expands (e.g. iDavel 
l2008h . 

The impact of halo quenching is clearly seen at masses 
> 3 X 10^° Mq. Below this transition mass, there is a mini- 
mal difference in the quenched and non-quenched relations, 
never exceeding a difference of more than 0.2 dex in the 
median: at these masses, star-forming galaxies still domi- 
nate the number densities of galaxies. Above this transition 
mass, the sSFR drops quickly towards s mall values, which i s 
qualitatively similar to observations (e.g. lSalim et al.ll2007^ . 
The transition mass evolves very little with redshift, being 
only marginally higher at z ~ 2, which reflects the fairly 
constant relation between halo mass and stellar mass in 
these simulations. At all redshifts, the highest mass galax- 
ies are still forming stars at a non-trivial rate, though much 
reduced relative to the no-quenching case. This partly re- 
flects the infall of satellites that have not fully quenched 
onto the central (quenched) object. Tracking the evolution 
of these objec ts in detail is a dif ficult problem particularly 
for SPH (e.g. lAgertz et al.|[2007l ). and hence the upturn in 
sSFR at the highest masses should be taken with some 
caution. Nevertheless, it is clear that quenching as imple- 
mented in our simulations has a dramatic impact on the 
sSFR of galaxies above the transition mass of 3 X 10^°Mq, 
where many galax y properties are observ ed to undergo a 
strong transition (|Kauffmann et al] |2004| ). In our model 
there are very few galaxies above this mass on the main 
star-forming sequence, and SAMs with similar quenching 



mech anisms suggest this be havior disagrees with observa- 
tions ISomerville et al.ll2008h 

Particularly noteworthy are the intermediate-mass 
galaxies that we argued earlier are "filling in" the gap of 
the red sequence at stellar masses of ^ 10^° Mq. As we will 
argue in § O these galaxies are a combination of large satel- 
lites and centrals that are below the quenching threshold but 
have been quenched owing to living in a hot large-scale envi- 
ronment. As these galaxies become more numerous towards 
lower redshifts, the transition between quenched and un- 
quenched is significantly blurred, such that many quenched 
galaxies are present at masses below the nominal transi- 
tion mass. Hence the argument that the observed gradual 
transition from unquenched to quenched galaxies rules out 
a model based on halo mass quenching (which our model 
approximates) may not be justified. While our simulations 
qualitatively resemble observations, it is worth noting that 
such environmental quenching processes are not as robustly 
modeled as we would like in our simulations, and hence a 
more quantitative comparison will have to await improved 
modeling of the interactions of galaxies moving through hot 
halo gas. 



4 GROWTH ALONG THE RED SEQUENCE 
4.1 Mergers since quenching 

We have seen that a galaxy does not stop growing and evolv- 
ing once it reaches the red sequence (see Figure [3)|. Its exist- 
ing stellar population ages, that population sheds mass via 
stellar winds, and it obtains new stars (and gas) via merg- 
ers. Here we show that the most massive galaxies, which are 
quenched at early times {z > 1), often grow by factors of 
a few between quenching and the present day. In contrast, 
low-mass red galaxies (mostly satellites) tend to lose mass 
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Figure 5. Growth of galaxies along the red sequence. In our model, the most massive galaxies typically grew by factors of a few since 
being quenched at 2 ~ 1, and that mass growth is dominated by minor mergers. Top Left: median (solid line), 10th and 90th percentiles 
(dashed) of the redshift at which red sequence galaxies were quenched (zquonch) in different mass bins. Bottom Left: fractional mass 
growth of red sequence galaxies since Zquonchi ^ ^ function of stellar mass at z = 0, among galaxies quenched before z = 0.5. Galaxy 
points are colour-coded by Zquenchi a-nd galaxies with higher Zqucnch have larger symbol sizes. Top Right: the fraction of galaxies per 
bin that have undergone at least one major merger since being quenched, among galaxies with Zquench > 0.5. The fraction for galaxies 
with 0.5 < Zquench < 1-0 (dotted) is typically lower than that for galaxies quenched earlier, z^^^-^^^ > 1.5 (dashed). Bottom Right: the 
mass-weighted average merger mass ratio, among galaxies with Zquen^h > 0-5- The mean merger is always a minor merger (< 1/3, thin 
dotted line). Galaxies quenched at high redshifts (dashed line) all have a merger mass ratio of 0.2, while those quenched at somewhat 
lower redshifts (dotted) show stronger variation with mass. 



through stellar evolution and stripping without undergoing 
significant mass growth via mergers. 

The top left panel of Figure [5] shows the distribution 
of the redshift at which galaxies are quenched, Zquonch, as 
a function of the final stellar mass at « = 0. We plot the 
50th (solid), 10th, and 90th percentiles (dashed) of the dis- 
tribution in each bin. The quenching redshift is defined as 
the redshift when a galaxy first appears red in the colour- 
magnitude diagram, with an evolving dividing line between 
red and blue galaxies as described in H3.ll Low-mass galaxies 
on the 2 = red sequence, which are predominantly satel- 
lites, are quenched at fairly late times (generally z < 1). 
Massive galaxies > IO^^A/q are quenched earlier, typically 
at 2 > 1 and beyond. 

The bottom left panel of Figure [S] shows, as a function 
of the 2 = stellar mass, the fractional mass growth of 
red sequence galaxies between the redshift of quenching and 



the present day, y = Mstollar(2 = 0)/Mstcllar(z = ^quenched)- 

Since many galaxies are quenched at late times, and such 
galaxies have not had time to undergo mergers, we only 
plot the mass growth for galaxies with Zqucnch > 0.5. This 
gives a long enough time baseline 5 Gyr, much longer 
than the halo dynamical times) to assess red galaxy growth 
via mergers. Galaxy points are colour-coded by ^quench, as 
indicated on the right side of the panel. 

Departures from the y = 1 line (i.e. no change in stel- 
lar mass, dotted line) increase with quenching redshift: The 
blue points that are the most recently quenched galaxies 
typically lie near y = 1, and in fact are slightly below owing 
to stellar mass loss. At higher quenching redshifts, galaxies 
span a wider range in y. In some cases, galaxies lose mass 
catastrophically by tidal interactions with other galaxies, 
and end up near y — 0. A significant proportion of galax- 
ies quenched at earlier times with M > IO^^'^A^q, how- 
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ever, have grown substantially, by factors up to ~ 3. Since 
these massive galaxies live in high-density environments, and 
they've been quenched for longer, they have had more oppor- 
tunities to acquire mass via mergers than their lower-mass 
counterparts. Hence y reflects the amount of mass growth 
by mergers in massive galaxies since the time of quenching, 
and in fact underestimates the growth since stellar mass loss 
also occurs. 

To further explore the nature of post-quenching mass 
growth in our model, we use our merger tree to separate 
major mergers (those with a mass ratio larger than 1:3) 
from minor mergers (any merger with a lower mass ratio). 
We include in the minor merger component accretion of ex- 
isting stars that live in poorly-resolved galaxies, but this has 
little impact on the results. The top right panel of Figure [S] 
shows the fraction of z = red sequence galaxies that have 
undergone at least one major merger since being quenched. 
Again, we restrict this analysis to galaxies quenched before 
z = 0.5. We also split these galaxies further by quenching 
redshift into high- (dashed, Zquonch > 1.5) and low-redshift 
(dotted, 0.5 < Zquonch < 1.0) bins. 

Figure [5] shows that post-quenching major mergers are 
generally rare, happening in typically only ~ 10% of the 
cases at most masses, with a mild peak up to ~ 20% at 
M* ~ 10^^ Mq. Galaxies quenched at the earliest times have 
major merger fractions as high as 30%. At higher masses we 
begin to be limited by our simulation volume, as such galax- 
ies become sufficiently rare for mergers between them may 
be underrepresented. Nevertheless, since many red sequence 
galaxies have grown by much more than 20% in mass (even 
while undergoing stellar mass loss), this result immediately 
suggests that major mergers play a sub-dominant role in the 
overall growth of massive red galaxies. 

With little or no in situ star formation, the growth of 
these galaxies must therefore owe predominantly to minor 
mergers. The bottom right panel of Figure [S] shows the mean 
mass- weighted merger mass ratio as a function of z = stel- 
lar mass. By mass-weighting, we can assess the contribu- 
tion that mergers of a given mass ratio make to the overall 
growth of passive galaxies. Again we only include galaxies 
that became red at z > 0.5, and we also split the sample 
into high- and \ow-Zq„ench sub-samples. The typical merger 
mass ratio for red sequence galaxies in our simulations is 
~ 0.2, below the commonly- used 1/3 cutoff for major merg- 
ers (thin dotted line). Our typical value agrees with that ob- 
tai ned in high e r-reso lution individual galaxy re-simulations 
by lOser et al.l (|201ll ). Galaxies quenched early on all have 
roughly the same mean merger mass ratio, whereas galaxies 
quenched at later times show variation with z — stellar 
mass such that the most massive and least massive galaxies 
have higher mass ratios than those in the middle. At all stel- 
lar masses, red galaxy mass growth is dominated by minor 
mergers. 

In summary, massive red sequence galaxies in our sim- 
ulation grow by factors of up to several in mass between 
the time they are quenched (typically z ~ 1) and redshift 
0. Lower mass passive galaxies (both satellites and centrals) 
tend to be quenched at later epochs, and grow less, typically 
even diminishing in mass owing to stellar evolution. Mass 
growth at all masses is dominated by minor mergers with a 
characteristic ratio of 1:5, with at most a 30% contribution 
from major mergers. 
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Figure 6. The ratio of size growth to mass growth for red se- 
quence galaxies, since 2 = 1, as a function of their stellar mass 
at 2 = 0. A solid Une indicates the median. Points are colour- 
coded by the local galaxy density. More massive galaxies, which 
live in denser environments, show more size growth than lower- 
mass ones. In some rare cases the galaxy size may evolve several 
times (~ 3) more than the mass. This is roughly in line with 
observations, where the ratio is typically ~ 2 for evolution since 
z ~ 2. 



4.2 Size evolution 

Observations indicate that massive passive galaxies are more 
compact at high redshifts than at z ~ by factors of 
up to 6 or more (iToft et al.ll2007l : Ivan der Wei et al.ll2008l : 
Ivan Dokkum et "aL 20081 '). While observational biases pre- 
sented an initial concern, recent work suggests that much 
of the evolution must be physical, and probably arises from 
a combination of several physical processes. These processes 
include adiabatic expansion associated with stellar mass 
loss, major mergers , minor mergers, and evolut i on in mass- 



to-light r atios (e.g 



[2009; Hop kins et al 



Bovlan-Kolchin et akllioo^ : iNaab et all 
" 2010t ). 



Our simulations make specific predictions about the dis- 
tribution of merger mass ratios that contribute to red galaxy 
growth. Therefore, even though the spatial resolution of our 
simulations is not well-suited to study the internal structure 
of individual galaxies, we can use simple analytic merger- 
based estimates to approximate the size growth that our 
simulated galaxies might undergo. 

We construct a simple model for t he effects o f merg - 
ers following the formalism presented in lNaab et al.l (|2009h . 
which arises directly from the virial theorem plus energy 
conservation. This is appropriate for a purely collisionless 
system, so we implicitly ignore any (presumably minor) ef- 
fects of dissipation. Given a stellar accretion event (i.e. a 
merger) of mass Ma onto a galaxy of mass Mg, the ratio of 
final to initial size of the galaxy is given by: 



LI 



(1) 



Here, rj — Ma/Mg is the fractional mass increase, and 
e =< «^ > / < Ug > is the ratio of the two galaxies' stellar 
populations' mean squ are velocity d ispers ions. Using high- 
resolution simulations, iNaab et al.l (|2009r ) has shown that 
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this formula accurately describes the size evolution of early- 
type galaxies undergoing hierarchical mergers. 

To estimate size growth, we must therefore estimate 
the mass ratio and the velocity dispersion ratio between 
the satellite (i.e. less massive) and central (more massive) 
galaxies. We know the ratio of stellar masses directly from 
our merger tree. For the velocity dispersion, since we can- 
not accurately resolve it in our simulations, we instead em- 
ploy an observationa l result relating it to stellar mass from 
iGallazzi et al] l|2006l ): logcr = -0.895 + 0.286 log M^teiiar- 
With this relation, we obtain e = {oajogf = [MalUgf-^"^ . 
Thus, for each merger event, we can calculate the fractional 
size growth of the remnant galaxy. Taking each z = red 
galaxy, we consider its entire merger history since being 
quenched. We account for the size growth of each of its 
merger events and multiply them to find the galaxy's to- 
tal fractional mass growth. 

We show the results of this exercise in Figure [6] which 
shows the size growth relative to mass growth (rfMi/riMf) 
as a function of stellar mass. We plot this quantitiy instead of 
merely size growth because massive red galaxies in our sim- 
ulations may not grow the correct amount of stellar mass, 
as hinted by the number density discrepancies in i]3.1l We 
account for this difficulty by comparing the size growth rel- 
ative to the mass growth of galaxies, which is a reflection of 
the prevalence of minor vs. major mergers. 

Based on this simple modeling, minor mergers appear 
to be numerous enough to drive the factors of several in size 
growth of passive galaxies over time. Galaxies with stellar 
masses > 10^^ Mq show typical size growth a factor of 2 
greater than their mass growth, indicating that a series of 
mi nor mergers tends to puff up these galaxies, as argued 
bv lNaab et al] (|2009D . Lower mass galaxies have not grown 
as much since moving onto the red sequence, since as we 
showed in Figure[5]they undergo less merging. These results, 
using a simple analytic model, are broadly consistent with 
the idea that minor mergers drive the rapid size growth seen 
in massive red galaxies since high redshift. 

In Figure[6]we also colour-code galaxies by their [z — 0) 
local galaxy densities (which we discuss further in the fol- 
lowing section). We note that all galaxies shown are in 
relatively dense regions compared to blue galaxies. Den- 
sity is correlated with stellar mass for these passive galax- 
ies ( H5.2|l . so that the galaxies in high-density regions (red 
points) lie on the right side of the plot. Since size growth 
(relative to mass growth) appears to increase with stellar 
mass in Figure |S1 galaxies in denser regions tend to un- 
dergo more size evolution than those in less-dense regions. 
At fixed mass ~ 10^^ Mq there is a weak correlation be- 
tween size growth and density. This qualitatively agrees 
with recent observations indicating that galaxies in high- 
density regions must undergo more evolution in the mass- 
size relation sin ce z = 1 than their counterparts in lower- 
density regions (|Raichoor et al.ll2011^ . Physically, based on 
our merger-driven size growth model, this suggests that red 
galaxies in the densest regions (i.e. galaxy groups in our 
simulations) tend to undergo a higher proportion of minor 
mergers. 

In summary, a simple model of red galaxy size growth 
driven by mergers shows that massive galaxies with > 
10" Mq in our simulation typically have grown in size by a 
factor of 2 relative to their mass growth since z — 1. The 



degree of size growth increases with stellar mass, and galax- 
ies in denser regions grow by larger factors. These results 
are broadly consistent with observations as well as previous 
simulation results. 



5 ENVIRONMENTAL FACTORS 

5.1 Bivariate dependence of environment on 
colour and A/stdiar 

We have shown that the massive and low-mass ends of the 
red sequence grow first in our model - not until z < 0.3 
do galaxies around IO^^Mq appear in comparable numbers 
along the red sequence. Surprisingly, some of these galax- 
ies are identified as central galaxies in a halo that is below 
the quenching threshold, rather than as satellites within a 
quenched halo. In this section, we examine how environmen- 
tal factors drive this behaviour. 

We illustrate these environmental effects in Figure[7l On 
the left we show a colour-mass diagram where the galaxies 
are colour-coded by their local galaxy density. The resulting 
trend is clear: the very highest and low-mass red galaxies 
live in the densest regions, while intermediate-mass passive 
galaxies (M« ~ 10^" — IO^^Mq) live in lower-density re- 
gions. This trend matc hes that seen in SPSS observations, 
at least qualitativelv (|Hpgg et al] l2003l : iKauffmann et ahl 
|2004 iBlanton et"ai] [loO^TPhysically, this trend emerges 
because massive galaxies with M* > 10^"'^ Mq, once mas- 
sive enough, can get to the red sequence "on their own" by 
forming their own hot halos even in (relatively) low-density 
environments. Massive galaxies that are quenched at early 
epochs live in the highest overdensities, and thus grow to be 
the most massive galaxies at z = by accreting other galax- 
ies (see Low-mass (< IO^^'^Mq) galaxies can typically 
only be quenched by living in the hot environment of more 
massive neighbors. 

In the right panel of Figure[7]we show only central galax- 
ies colour-coded by the distance to the nearest halo with 
a mass > IO^^Mq. This is roughly the halo mass where 
hot halos will form independently in our simulations. Gray 
crosses are for galaxies within 10h~^ kpc of the centers of 
massive halos, and coloured points are for central galaxies 
in lower-mass halos. Note that there are many such galaxies 
along the red sequence at intermediate masses. We colour- 
code each point based on the ratio of its distance from the 
center of the nearest massive halo to that halo's virial radius. 

The first clear difference when looking only at centrals is 
that there is a deficit of A/, ^ 10^" red central galaxies, 
which as we argued in >i 3. 3 1 arises because these galaxies are 
predominantly quenched satellites that have grown mostly 
at late times. 

However, there are still a fair number of central galax- 
ies that are along the red sequence at intermediate to low 
masses. Such galaxies at low mass tend to be fairly close 
to a quenched halo (i.e. the points are blue), and at least 
some of the intermediate mass ones as well. This suggests 
a larger environmental effect - these galaxies are part of 
super-group structures around higher-mass galaxies, whose 
hot (albeit unvirialised) gas can be felt beyond the virial ra- 
dius of member galaxies. Alternatively, these galaxies may 
have passed within _Rvir of a hot halo on an eccentric orbit 
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Figure 7. Left CMD (at 2 = 0) with galaxies colour-coded by the local galaxy density, which is computed by counting the number 
of galaxies within a radius of 1 Mpc. The high- and low-mass ends of the red sequence live in dense environments, whereas ~ L* red 
sequence galaxies tend to live at lower densities. Right Central galaxies are colour-coded by the distance to the closest massive halo 
(right). Massive halos are defined by a mass of 10^^ Mq, as found using a spherical overdensity algorithm. This is roughly the mass at 
which a halo is dominated by hot gas (> 250, 000 K). Grey crosses are galaxies at the centers of their (massive) halos, and blue and red 
points are outside the virial radius of the nearest massive halo. Some passive galaxies < IO^'^'^Mq are quenched despite living outside 
-Rviri suggesting the influence of hot gas in unrelaxed super-group environments. Quenched galaxies at > 4i?vir live in halos just below 
where scatter in the Mij^iq— hot gas relation allows them to be hot-gas-dominated. 
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Figure 8. Contours of red fraction in stellar mass-overdensity 
space in our quenching simulation. Galaxies are increasingly red 
at both higher masses and overdensities. Small galaxies at high 
overdensity are red because they live within or close to larger 
halos dominated by hot gas, while larger galaxies are red because 
they are typically centrals living with hot halos maintained hot 
via our quenching feedback mechanism. 

and been quenched during close passage. This result suggests 
that environmental factors beyond the halo mass may play 
some role in quenching of star formation. Observations have 



suggested a role fo r environment (beyon d halo mass) in some 
galaxy properties l|Cooper et al.ll2008al lbl. I2OI0I ) , though it is 
typically secondary to galaxy mass. 

Figure [8] shows another representation of the relation- 
ship between red fraction and environment. Here we plot 
contours of red fraction (in colors from blue to red) in the 
plane of stellar mass versus local galaxy overdensity (de- 
termined u sing Smooth). Th is figure can be compared to 
Figure 6 of iPeng et all (|2010| ). who argued for two different 
quenching mechanisms, namely "environment quenching" of 
all galaxies above a given overdensity, and "mass quenching" 
above a given stellar mass. 

Our simulation quantitatively reproduces the observed 
trends almost exactly, a lthough it has far fewer galaxies 
than in the SDSS data of iPeng et el\ (I2OI0I ). and hence the 
trends are less smooth and there are regions in this space 
where this simulation yields no galaxies (black). Broadly, 
our models reproduce the strong increase in red fraction to 
both high mass and high overd e nsity, and the "boxy" con- 
tour shape noted by IPeng et ah! (|2010l ) is indicative of these 
quenching mechanisms being mostly separable. In our simu- 
lations, environment quenching is predominantly associated 
with satellites at lower masses living in hot halos, although 
as discussed above there are a fair number of small cen- 
trals living in dense environments that are also quenched. 
Meanwhile, mass quenching is well reproduced by our hot 
gas threshold for quenching, which translates roughly into a 
halo mass threshold. 

In summar y, our simulations are consistent with the in- 
terpretation bv lPeng et al.l (|2010l ) that there are two sepa- 
rable quenching mechanisms that deliver galaxies to the red 
sequence, associated with environment and mass. But our 





Figure 9. Top row: Mass-density relation for all (solid), blue (dotted), and red (dashed) galaxies at 2 = 0, 1, and 2 (increasing left-to- 
right). On average (in the "all" sample), galaxies below 10^"'^ Mq consistently live in low density environments, with a sharp rise in local 
density at higher masses. Red galaxies at low masses live in high-density environments, producing a pronounced dip in the red galaxy 
mass-density relation at ~ IO^^'^Mq. Redshift evolution to 2 = 1 is negligible, and weak at higher z. Bottom row: Colour-density 
relation for all (solid), low-mass (dotted), and high- mass (dashed) galaxies at three redshifts. On average, blue galaxies live in low-density 
environments and red galaxies in high-density environments. As expected from the mass-density relation in the top row, massive galaxies 
tend to live in denser environments than low-mass galaxies, except perhaps for massive red galaxies. While the overall colour-evolution 
of galaxies is apparent in the different panels, the basic trends appear roughly the same at all z. 



simulations further suggest that both mechanisms are ulti- 
mately the result of red galaxies living in regions dominated 
by hot gas, which occurs at both high overdensities and high 
masses. 



5.2 Evolution of the mass-density and 
colour-density relations 

We show the evolution of the stellar mass-density relation 
and the colour-density relation in Figure [S] Here the local 
galaxy density is found using Smooth. 

The mass-density relation (top row of Figure ^ for all 
galaxies shows a very shallow rise in density with mass for 
stellar masses below 10^°-'^ Mq. Ah ove this mass, density 
rises steeply. This behavior agrees qualitatively with obser- 
vations at 2 = 0, in which a fiat slope at lo w luminosities be- 
comes very steep for the brightest galaxies (|Hogg et al.ll200 j : 
iBlanton et ahlboOSl ). The trends for red and blue galaxies, 
however, are distinct: Blue galaxies (dotted lines) tend to 
live in low-density regions and follow the trend for all galax- 
ies at the low-mass end (since they dominate by number). 
Low-mass and high-mass passive galaxies live in dense re- 
gions, while passive galaxies with M« 10^" — 10^^ Mq live 
in lower-density regions, as discussed in § 15.11 

For these trends, there is negligible evolution out to 



z = 1 (top middle panel). Measurements of the z — 1 mass- 
and colour-density relations have generated some debate 
in the literature, with recent in dications that the colour- 
density rel ation persists to z ~ 1 | Cooper et al."2006', '2007|; 
lElbaz et al . 2007; Scodc ggio et al.i2 009: Cooper ct al. 201(J), 
in accordance with the predicted lack of evolution in Figure 
[9l At z = 2 (right panel), however, massive blue galaxies 
live in higher density environments than massive red galax- 
ies. This suggests that the colour-density relation may be 
inverted for the most massive galaxies at high redshifts in 
our models. There are two cautionary notes for this, how- 
ever: first, this result may be sensitive to the division be- 
tween red and blue galaxies: Figure [1] shows that many blue 
galaxies are moving toward the red sequence in the z — 2 
panel. Second, our quenching criterion may not be capturing 
the correct physics at these redshifts owing to cold streams 
penetrating hot halos. Therefore we do not view this high- 
redshift colour-density inversion as a strong prediction of 
our models at this time. 

The colour-density relation (bottom row of Figure [5| 
shows similar behavior to the mass-density relation: den- 
sity is constant for blue colours (low values of g — r), then 
rises steeply for red colours (higher values of g — r). When 
split into high (dashed) and low (dotted) mass galaxies (at 
M, = 1O^°-^M0), we see that low-mass galaxies follow the 
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same trend as for all galaxies. High-mass galaxies, on the 
other hand, live in higher densities at the blue end, and 
lower densities at the red end. There are only a handful of 
blue galaxies with M > 10^°'^ Mq, and as expected these 
tend to live in higher density regions than less massive blue 
galaxies. The fact that the red end of the high-mass curve 
at 2 = is at lower densities than most passive galaxies 
(solid line) is another manifestation of the two-peak shape 
in the mass-density relation for the red seqence (Figure [T)): 
Red galaxies above 10^°'^A/q are a mixture of very mas- 
sive galaxies in high-density regions and ~ AI* galaxies that 
are relatively isolated. In contrast, low-mass passive galax- 
ies uniformly reside in high-density regions, since that's the 
only circumstance under which they become red. 

The overall colour-density relation at z = 1 (bottom 
middle panel) is not significantly different from that at 
z — 0. The relation does show some minor evolution for 
massive galaxies, however. At the red end at z = 1, high- 
mass galaxies lie on the overall trend, unlike the z — case 
where massive red galaxies are in less dense regions. This 
evolution suggests that many massive galaxies that move 
onto the red sequence between z = 1 and z = are rel- 
atively isolated, and reflects the increasing contribution of 
intermediate-mass red galaxies with time. 

In summary, although our quenching model is most di- 
rectly tied to hot gas fraction and thus halo mass, envi- 
ronment plays a correlated role in quenching galaxies in 
our simulations. Since our quenching is directly linked to 
the proportion of hot ga s in a halo and thus halo mass 
jBirnboim fc Dekej l2003l : iKeres et all l2005l : iGabor et all 
and halo mass is correlated with environment in a 
CDM universe, this naturally leads to galaxies in the dens- 
est environments being quenched first. Environment plays a 
key role for a significant fraction of intermediate- mass galax- 
ies that are not in quenched halos themselves but live near 
quenched halos; these become quenched owing to living in 
an overdense region pervaded by hot gas beyond the virial 
radius of massive halos. The latest-arriving members of the 
red sequence near M* typically live in less dense regions, 
namely isolated halos of M > IO^^Mq. Evolut ion in the 
mass-density and colour-density relations is weak in this 
model and broadly consistent with data, perhaps showing 
an inversion in the colour-density relation only at the high- 
est redshifts considered here. 



6 SUMMARY AND CONCLUSION 

We have analysed a cosmological hydrodynamic simulation 
to explore the evolution of the red sequence. A key feature of 
our simulation is a simple quenching prescription tied to the 
dominance of hot gas within the virial radius of a halo. In our 
quenching prescription, galaxies whose halos have > 60% of 
their baryons in hot gas are starved of star-forming fuel by 
continual heating of the halo gas. In a previous paper we 
showed that this simple model yields a distinct red sequence 
of galaxies whose number den sities generally agree with ob- 
servations (jCabor et al.ll201l] '). This model thus provides a 
tool for studying the emergence and growth of the red se- 
quence over cosmic time. While some details of our model 
are overly simplistic and in some cases are in confiict with 
observations, it serves as representative of a general class of 



quenching models where cessation of star formation is driven 
by the presence of hot coronae or is tied to halo mass. 

In our simulations, the first galaxies to turn red at 
2 > 2 are among the most massive in the universe at 
that epoch, with A/stdiar ~ 10^"'^~"Mq. Satellite galax- 
ies at low masses are concurrently starved of fuel by the 
same hot halo that triggers our quenching mechanism, and 
therefore low mass galaxies also appear on the red sequence 
at early times. Meanwhile, intermediate-mass galaxies at 
~ 10^'' — W^'^'^Mq remain predominantly blue until late 
epochs, producing a pronounced dip in the red galaxy mass 
function at high redshifts. This dearth of low- mass red galax- 
ies is consistent with observations at moderate (< 1) and 
high redshift (z > 2 — 3). At later times, the number of 
moderate-mass red galaxies catches up to the number of 
more massive ones, so that the dip in the red galaxy mass 
function is less pronounced by 2 = 0. 

Essentially all simulated passive galaxies with stellar 
mass > lO^^A^Q form from smaller passive galaxies that 
grow mass via mergers once on the red sequence. This mas- 
sive end of the red sequence does not emerge until 2 ~ 1.5, 
and it grows in number by factors of several at late times. 

Most of the mass growth of individual passive galaxies 
is due to minor mergers, with major mergers being rare and 
a sub-dominant component of mass growth. Based on a sim- 
ple analytic model, these numerous minor mergers can lead 
to size growth that is a factor of 2 or more times that of a 
galaxy's corresponding mass growth. This suggests that mi- 
nor mergers dominate the dramatic size ev olution of massive 
galaxies, in accorda nce with earlier work (jNaab et al ] |2009l : 
iHopkins et allbOlOl '). 

The simple quenching model explored here qualitatively 
reproduces observed trends among colour, stellar mass, and 
environment at 2 = 0. In detail, it also reproduces the ob- 
served mass-density trend for the red sequence: low-mass 
and high-mass passive galaxies live in the densest regions, 
whereas galaxies in-between (around M*) live in less dense 
regions. The mass-density and colour-density trends persist 
to 2 = 1, but the colour-density trend seems to invert for 
the most massive galaxies at 2 = 2. Quenching occurs at 
both high density, roughly independent of stellar mass, and 
high stellar mass, roughly independ ent of density. As a rgued 
from corresponding observations bv lPeng et al.l (|2010l ). this 
indicates two quenching mechanisms associated with envi- 
ronment and mass, which in our model both arise when the 
local environment is dominated by hot gas. 

Our analysis explictly avoids any considerations of 
galaxy morphologies. A key question that we therefore do 
not address is, why are most red and dead galaxies ellip- 
tical? While the conventional thinking is that it owes to 
major mergers th at can transform star-form ing spirals into 
elliptical galaxies ijMihos &: Hernauistlll99^ ). the simple ces- 
sation of star formation at sufficiently early epochs typically 
results in substantial later growth via dry (i.e. purely or al- 
most purely stellar) mergers. In this scenario, galaxies would 
quench first when entering a h ot halo (or its imme diate envi- 
ronment), become red spirals (|Bundv et al.ll2010h . and then 
undergo subsequent dry mergers that alter its morphology. 
Our models indicate that larger galaxies living in denser en- 
virons undergo more such growth, and so would be expected 
to have dynamical signatures indicating more dry merging 
such as boxier isophotes; these trends are qualitatively con- 
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sistent with observations. Whether the quantitative predic- 
tions of this model are in accord with detailed observations 
remains to be seen. 

More broadly, it appears that a quenching mechanism 
based on a simple hot gas fraction threshold (which is well- 
correlated with a halo mass threshold) is capable of repro- 
ducing a variety of observations not only at the present 
epoch but also out to ~ 1 and beyond. However, there 
are significant failures, such as a dearth of the most mas- 
sive galaxies, and a too-sharp cutoff in the blue galaxy mass 
function suggesting that both today and at high redshifts 
in the real Universe, some galaxies are able to be star- 
forming despite living in hot gas-dominated halos. Hence 
this model represents only a step towards understanding 
the physical processes driving red galaxy evolution. Compar- 
isons to observations particularly probing the intermediate- 
mass regime out to high redshifts can place strong con- 
straints on models such as these, and highlight additional 
physical processes that may need to be considered. 
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